Spatially nonuniform phases in the one-dimensional SU(n) Hubbard model 

for commensurate fillings 
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The one-dimensional repulsive SU(n) Hubbard model is investigated analytically by bosonization 
approach and numerically using the density-matrix renormalization-group (DMRG) method for 
n = 3, 4, and 5 for commensurate fillings f — p/q where p and q are relatively prime. It is shown 
that the behavior of the system is drastically different depending on whether q > n, q — n, or q < n. 
When q > n, the umklapp processes are irrelevant, the model is equivalent to an n-component 
Luttinger liquid with central charge c — n. When q = n, the charge and spin modes are decoupled, 
the umklapp processes open a charge gap for finite U > 0, whereas the spin modes remain gapless 
and the central charge c — n — 1. The translational symmetry is not broken in the ground state 
for any n. On the other hand, when q < n, the charge and spin modes are coupled, the umklapp 
processes open gaps in all excitation branches, and a spatially nonuniform ground state develops. 
Bond-ordered dimerized, trimerized or tetramerized phases are found depending on the filling. 

PACS numbers: 71.10.Fd 
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I. INTRODUCTION 

Recently, the SU(n)-symmetric generalization of the 
standard SU(2) Hubbard model! has been intensively 
studied theoretically^i^iii&i&Si Apart from its theoret- 
ical interest this model may mimic strongly correlated 
electron systems where the orbital degrees of freedom of 
d and / electrons play important role and these extra 
degrees of freedom are taken into account by consider- 
ing n-component fermions. On the other hand, ultracold 
gases in optical lattices may also be simulated by such 
multi-component models. 

The Hamiltonian of the model is usually written in the 
form 

N n 



U 



i=l (7=1 

N n 



1=1 (7,cr'=l 



(1) 



where N is the number of sites in the chain. The oper- 
ator c| (Cjg.) creates (annihilates) an electron at site i 
with spin a, where the spin index is allowed to take n 
different values. n i iy is the particle-number operator, t is 
the hopping integral between nearest-neighbor sites, and 
U is the strength of the on-site Coulomb repulsion. In 
what follows t will be taken as the unit of energy. 

The model behaves as an n-component Tomonaga- 
Luttinger liquid at generic fillings. Other type of behav- 
ior may appear at commensurate fillings due to umklapp 
processes. The possible phases, their nature and the crit- 
ical coupling where they appear have been studied in 
detail for two special commensurate fillings of the band, 
namely for half filling and 1 /n filling It is well 

established by now that the ground state is a fully gapped 
bond-ordered dimerized state in the half-filled case for 



any n > 2. Contrary to this, the ground state remains 
translationally invariant in the 1/n-filled case, and only 
the charge mode acquires a gap for U > U c . While As- 
saraf et al 3 argued that U c is finite, our recent numerical 
work-2 has suggested a much less, perhaps U c — crit- 
ical value above which multiparticle umklapp processes 
become relevant. 

It is worth mentioning that the SU(n) Hubbard model 
has a rich phase diagram in the attractive case, too. 10 : 11 
The one-third-fillcd SU(3) model has two distinct phases 
in the high-dimensional limit. In one of them the fer- 
mions form trions while in the other phase a color super- 
fluid state emerges. The existence of these phases is not 
yet settled in one dimension. 

In the present paper, the role of multiparticle umklapp 
processes will be further analyzed for general commensu- 
rate fillings / = p/q, where p and q are relatively prime. 
We try to establish under what conditions the umklapp 
processes can generate gaps in the charge or spin sectors, 
and when and how the translational symmetry is broken. 
To this end, partly analytic, partly numerical procedures 
will be applied. We will generalize the method used in 
Ref. to the one-third-filled SU(n) model to show an- 
alytically that the ground state cannot be spatially uni- 
form. It is trimerized at least in the large-n limit. The 
numerical work will show that in fact this trimerized state 
with gapped excitations exists for n > 3 already. 

In the numerical part, the length-dependence of the en- 
tropy of finite blocks of a long chain is studied. Recently, 
it has been shown that quantum phase transitions can 
be conveniently studied by calculating some measure of 
entanglement This can either be a 

local quantity, e.g., the concurrence^!, a global quantity, 
e.g., the fidelity, 22 or the entropy of a block of several 
sitesi 2 ^ As has been demonstrated recently^ the oscil- 
latory behavior of the block entropy can reveal the posi- 
tion of soft modes in the excitation spectrum of critical 
systems or the spatial inhomogeneity of gapped models. 
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This will allow us to demonstrate that at commensurate 
fillings / = p/q the type of ground state of the one- 
dimensional SU(n) models depends on whether q = n, 
q < n, or q > n. 

The paper is organized as follows. The oscillatory be- 
havior of the block entropy, the corresponding peaks in 
its Fourier spectrum, and their relationship to the known 
properties of the half-filled and 1/n-filled models are re- 
called in Sec. II, where some new results necessitating 
further studies are also given. An analytical investigation 
of the role of umklapp processes at commensurate fillings 
is presented in Sec. Ill and the possibility of spatial inho- 
mogeneity of the ground state is discussed. The numer- 
ical results for various fillings are presented in Sec. IV. 
Finally our findings and conclusions are summarized in 
Sec. V. 



II. OSCILLATORY LENGTH DEPENDENCE OF 
THE BLOCK ENTROPY 



If a finite block of length / of a long chain of N sites is 
considered, it is in a mixed state, even if the long chain 
is in its ground state. The mixed state can be described 
by a density matrix pn(1) and the corresponding von 
Neumann entropy is 

s N (l) = -Tr[p N (l)lnp N (l)]. (2) 

It is well know n 23 ' 25 that this entropy as a function of the 
block size grows logarithmically if the system is critical 
and the spectrum is gapless. In addition, the central 
charge c can be derive d 26 ' 27 from the initial slope of the 
length dependence of sn(1), 



Sn(1) = r In 
6 
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(3) 



where g is a shift due to the open boundary. It contains a 
constant term which depends on the ground-state degen- 
eracy and an alternating term decaying with a power of 
the distance from the boundary 28 ! 29 On the other hand 
for noncritical, gapped models, sjv(0 saturates to a finite 
value when I is far from the boundaries. 

Recently it has been pointed out by some of us 2 ^ that 
a wider variety of behavior may be found for the length- 
dependence of the block entropy. Namely, we have shown 
that in some cases oscillations may appear in sat(Z). This 
can be best analyzed by considering the Fourier trans- 
form 



~s(k) 



1 N 

-Y 



-ikl 



S N (l) 



(4) 



1=0 



for discrete wave numbers, k = 2tt j/N lying in the range 
(— 7T,7r). Since sn(1) — sn(N — I), s(k) is real. It has 
a large peak at k = 0, all other Fourier components are 
negative. Peaks in \s(k)\ carry information about the po- 
sition of soft modes or the spatial inhomogeneity of the 



ground state. More precisely if the amplitude of a peak 
at a nonzero wave number k* remains finite in the ther- 
modynamic limit, this indicates a periodic spatial mod- 
ulation of the ground state with wavelength A = 2n/k* . 
On the other hand, if a marked peak appears in \s(k)\ but 
its amplitude vanishes as N — > oo, this allows to identify 
the wave vector of soft modes in critical models. 

In a recent work 9 - we have shown that such oscillations 
appear in Sn(1) for the SU(n) Hubbard model as well. 
The periodicity depends on both n and the band filling 
/ = p/q. This is shown for the n — 3 and n = 4 models 
for the 1 / n-filled and half-filled cases in Fig. [1] for a large 
value of U (U = 10). 
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FIG. 1: Block entropy sjv(Z) of finite chains with N — 18 and 
N — 16 site, respectively, for n = 3 and 4 at fillings / = 1/n 
and / = 1/2 for U = 10. The solid line is our fit using Eq. (J3]|. 

In the 1/n-filled cases, sjv(7) increases logarithmically 
with the block length (and then goes down as I ap- 
proaches N). When every third or fourth values are 
taken, depending on the periodicity, these selected val- 
ues can be fitted to ((3]) as shown by the solid lines in 
panels (a) and (c). This indicates gapless behavior and 
gives c = n — 1. This is in agreement with the theoretical 
expectation, since the charge mode becomes gapped due 
to multiparticle umklapp processes and only the n — 1 
spin modes are gapless. A distinct behavior is found in 
half- filled systems, as seen in panels (b) and (d). The 
quantity sjv(Z) oscillates with period two, and if only ev- 
ery second point is taken, it seems to saturate beyond 
some block length, before decreasing again, indicating 
that the corresponding models are fully gapped. 

The finite-size dependence of the peaks of \s(k)\ ap- 
pearing at k* — 2kp — 2irf characterizing the oscilla- 
tion is shown in Fig. [2] It is seen that in the 1 /n- filled 
case the Fourier components at k* = 2ir/n vanish in the 
thermodynamic limit, while a finite value is obtained at 
k* = 7r for half-filled models. This corroborates our find- 
ing that the 1/n-filled SU(n) models are critical with 
a spatially uniform ground state while a gapped bond- 
ordered dimerized phase appears at half filling. 

We have done similar calculations for more general 
commensurate fillings of the band. Figure [3] shows the 
results obtained for n = 3, / = 2/5 as well as for n = 4, 
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/ = 1/3. When every fifth points are taken for the 
/ = 2/5 filled SU(3) model, they can be fitted to © 
yielding c = 3 for the central charge. This indicates that 
all modes, including the charge mode, are gapless. 
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FIG. 3: Same as Fig. [T] but for (a) n = 3, / = 2/5, N = 20 
and (b) n = 4, / = 1/3, N = 24. 



Such a fit does not work for the one-third-filled SU(4) 
model. To better see the difference \s(k)\ is consid- 
ered again. The amplitude of the Fourier component at 
k* = 47r/5, also displayed in Fig. [2] for the n = 3 model, 
vanishes in the N — * oo limit. On the other hand \s(k)\ 
remains finite at k* = 2ir/3 in the one-third-filled n = 4 
model. 

When the same calculations are repeated for the n = 5 
model at / = 1/2,1/3,1/4, and 1/5, peaks appear in 
|s(fc)| at k* = it, 27r/3, 7r/2, and 27r/5, respectively. As 
is seen in Fig. [H the amplitude of the peaks remains 
finite even when N — > oo in the first three cases, while it 
vanishes in the last case. 

These results indicate that the role of umklapp pro- 
cesses depends on the relationship between the number 
of components n and the relative primes p and q char- 
acterizing the commensurate filling. In what follows this 
problem will be studied first analytically in a bosoniza- 
tion approach and large-n expansion technique, and then 
numerically using the DMRG method. 



III. ANALYTICAL CONSIDERATIONS 

A. The role of umklapp processes: a bosonization 
approach 

Following the usual procedure we write the Hamil- 
tonian JTJ) in momentum space and linearize the free- 
particle spectrum around the two Fermi points (±fcp). 
The underlying assumption is that the low-lying excita- 
tions determine the physics of the system. Depending on 
whether the momentum of the fermions is close to +k-p 
or — &f, one can distinguish left- and right-moving parti- 
cles, and the interaction processes also can be classified 
on the basis of whether the incoming and scattered parti- 
cles are right or left movers and the momentum transfer 
is small (forward scattering) or large, of the order of 2kp 
(backward scattering). In a generic model the strength 
of the various scattering processes may be different. For 
the sake of simplicity we neglect chiral processes in which 
both particles move in the same direction before and af- 
ter the interaction, since they lead to the renormalization 
of the Fermi velocity only. 

One can recognize that at generic fillings, where um- 
klapp processes do not play a role, the forward and back- 
ward scattering processes can be interpreted as current- 
current interactions and their contribution to the Hamil- 
tonian density can be conveniently rewritten using Dirac 
fermions- in the following short form (automatic sum- 
mation for the repeated indices is understood): 

Hint(x) = \g ailJ2IJ3ai 4! IJl {x)-i ll 'ip l y 2 {x)4) l y 3 {x)^ ll ^ ni {x). 

Here Oi denote the spin indices that can take the val- 
ues l,...,n, 7^ with fx = 1,2 are the Dirac matrices, 
in our case the standard Pauli matrices (a x , <r y ), and 
V>(x) = ip^(x)^i. While the Hubbard model contains a 
single interaction parameter U, the couplings ga^a^ai 
may be different for physically different processes in more 
realistic models. In the renormalization-group treatment 
we will assume this to be the case. It is assumed, how- 
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ever, that the spin of the fermions does not change in the 
scattering process and the couplings are symmetric un- 
der the exchange (a±, 02) «-> ((T3, (T4). If the fermion field 
ij) a (x) is decomposed into left- and right- moving compo- 
nents according to 



_ (R a {x) 



$„(x) = (4(a0, , (6) 



the usual backward- and forward-scattering processes are 
in fact recovered. In the standard g-ology^L notation 
g aa > a i a is denoted by -g x , and g aaa i a > by g 2 . 

The well-known renormalization-group (RG) equa- 
tions, the (3 function can be written for these scattering 
processes in a short form^ 

9 In g ai 020304 _ a _ 

9 In A' /A ~~ "o-l 0-20-3 0-4 — 9cr 1 a i cr 3 ij j gcr i a2cr j a i . 

9<j\ai<ijai9<Jia2cr%tTji 

where A is the cut-off parameter. These RG equations 
have been analyzed earlier— t 7 and it was found that the 
backward-scattering processes scale out at generic fill- 
ings in the SU(n) Hubbard model and for this reason 
this model is equivalent to an n-component Luttinger liq- 
uid in this case. The Hamiltonian can be diagonalized 8 
and the excitation spectrum can be determined exactly 
in bosonic phase-field representation. — There is one sym- 
metric combination of the phase fields with different spin 
indices, this is the so-called charge mode: 



1 - 

-7=22<t>cr( x )> 



(8) 



while the n—1 antisymmetric combinations give the spin 
modes: 



^ms{x) = 



1 



yjm[m + 1) 



} J <t>a{x) - TO0 m+ l(x) 



(9) 



with m = 1, . . . ,n — 1. 



Similarly to the spin-charge separation in the two- 
component Luttinger model, one finds complete mode 
separation. The Hamiltonian density is the sum of the 
contributions of the individual modes, 



(10) 



where j = c, Is, 2s, . . . , (n — l)s. Each term has the usual 
bosonic form: 

Bj(x) = ^ fa [Tl^x)} 2 + ld x Mx)f) , (11) 

where n,-(x) is the momentum canonically conjugated 
to 4>j{x). The renormalized velocities and the Luttinger 



parameters can be given in terms of the new couplings 
g2-j appearing after diagonalization^ in the spin indices: 

(12a) 
(12b) 



In a finite system, where the momentum is quantized 
in units of 2ir/L, the excitation spectrum of the Luttinger 
model can be written as: 




E = huj — f 'ir j 

j 



Ai + A 3 _ 



(13) 




where n± are integers describing the bosonic excitations 
and 



(14) 



where SNj is the change in the number of particles in the 
jth channel, and similarly Jj describes the current in the 
jth channel created by adding particles to or removing 
them from the branches of the dispersion relation. 
The total momentum is given by 



P = hk F J c + h— («+ - n°_ + A 3 + - A{ 



L 



(15) 



Thus soft modes appear not only at zero momentum but 
at integer multiples of 2kp, too, since the charge current 
J c is an even number if the total charge is conserved. 
Since 



2fa 



Nl2w 
n L 



(16) 



where N® is the number of particles in the system, and 
the filling of the band is / = N^/nN , the position of 
these soft modes depends on the filling only, 2fcp = 2%f. 

We know that the usual umklapp processes, scattering 
of two right movers into left-moving states or vice versa, 
that were neglected so far, are relevant in a half-filled 
system for any n££ Multiparticle umklapp processes may 
become relevant at other commensurate fillings. To see 
what kind of processes are allowed we have to take into 
account that the total quasimomentum transferred in an 
umklapp process has to be an integer multiple of 27r. If 
the band is / = p/q filled and consequently fcp = p^/q, q 
particles have to be scattered from one Fermi point to the 
opposite one to satisfy this condition. This is illustrated 
in Fig. [5j Since there arc n different types of fermion in 
the SU(n) model, and the Hubbard interaction is local, 
only such umklapp processes are allowed by the Pauli 
principle in which the number of scattered particles is less 
than or equal to n. Thus q-particle umklapp processes 
are forbidden in our model when q > n. 



5 



n = 3 
k F = tt/2 




gapless, thus c — n for {/ > 0. The expected behavior 
for different cases is summarized in Table |U 



n = 3 
= 7r/3 




FIG. 5: Umklapp processes in the SU(3) symmetric Hubbard 
model, (a) Any two of the three types of fermion can be 
scattered in the half-filled case, (b) Fermions with the three 
possible spin orientations participate in the scattering process 
in the one-third-filled case. 



The role of multiparticle umklapp processes in the con- 
ductivity of the SU(2) Hubbard model has been stud- 
ied in Ref. 33]. A different aspect, whether the charge 
and spin modes are coupled or not by the possible mul- 
tiparticle umklapp processes has been considered for the 
SU(n) models The g-particle umklapp processes can be 
described in the bosonic phase-field representation by 



H\j(x) = .93 



dx^2 cos [ 2 (0<ri (x)+... + (j} aq (x)) 
J {*«}' 



(17) 

The phase fields appearing here are the phases of the 
bosonic representation of the particles participating in 
the scattering process, and {cr,}' indicates that all spin 
indices have to be different. 

This Hamiltonian can be expressed in terms of the 
phase fields corresponding to the charge and spin modes, 
making use of the inverse of ([5]) and ©■ It is clear that 
if q = n, only the symmetric combination of the phase 
fields, that is the charge mode defined in (j8j) appears in 
(j!7[) . This means that in this case the g-particle umklapp 
processes couple to the charge mode only. The Hamilto- 
nian density of the charge mode is identical to that of the 
well-known sine-Gordon model which has a fully gapped 
excitation spectrum. Thus the charge mode becomes 
gapped for finite U. The n — 1 spin modes are not in- 
fluenced by the umklapp processes in the 1/n- filled case. 
They remain gapless and the central charge is c = n — 1. 

When q < n, the sum of the q phase fields correspond- 
ing to the q particles required for umklapp processes, will 
contain various combinations of the n boson fields, lead- 
ing to a mixing (coupling) of the charge and spin modes, 
thus opening gaps in all modes. The model becomes non- 
critical for U > U c . 

When on the other hand q > n, umklapp processes 
are forbidden by Pauli's exclusion principle when the in- 
teraction is local. The charge and spin modes remain 





n 


c phase 


k* 


q = n 


any n 


n - 1 C0S(ri - 1) 


2irp/n 


q < n 


n=£2 


coso 


2-Kp/q 


q > n 


any n 


n ClS(n-l) 


2np/q 



TABLE I: Central charge c and the type of phase charac- 
terized by the number of soft modes in the charge and spin 
sectors (CxSy) for the p/g-filled SU(n) Hubbard model. 



B. Spatial inhomogeneity in the large-n limit of the 
SU(ti) Hubbard-Heisenberg model at one-third 
filling 

We have seen that the spectrum is fully gapped when 
q < n. The question naturally arises whether the opening 
of a gap at multiples of k* = 2kp is related to a breaking 
of the translational symmetry, an instability against the 
formation of a spatially inhomogeneous state with the 
corresponding wave number. 

To analyze the stability of the homogeneous state we 
generalize the procedure used by Marston and Affleck^ 
to the one-third-filled case in the large-C/ limit. When n 
is an integer multiple of 3, the number of fermions sit- 
ting on each site is an integer in a homogeneous sample, 
and a finite energy is needed to add an extra particle. 
This energy gap at kp may imply a tripling of the spatial 
period in the ground state. To search for this spatial in- 
homogeneity a more general model, the SU(n) symmetric 
generalization of the Hubbard-Heisenberg model will be 
considerd. Its Hamiltonian is 



A 



^ = E 



j 



E (4^+1,^(4+1,^^,) (is) 



2-i 



CT=1 ^<T=1 ' 



The chemical potential is shifted to zero at one-third fill- 
ing, and the Hubbard coupling U and Heisenberg cou- 
pling J are rescaled by 2 /n so that the spacing of the en- 
ergy levels remain the same as n increases. We note that 
the usual Jj2<ij> ^i^j nearest-neighbor Heisenberg in- 
teraction breaks up into three terms in fermionic repre- 
sentation of the spin operators. One of them only shifts 
the chemical potential, another corresponds to nearest- 
neighbor Coulomb interaction (that is unimportant in the 
large U limit), therefore they are neglected. 

The equilibrium state will be determined from the 
minimum of the free energy that can be derived from 
the partition function of the system. In functional inte- 
gral formalism the partition function can be expressed 
with the Lagrangian of the model. At finite tem- 
peratures the imaginary time Lagrangian is L[c, c^] = 



6 



£i a ( c l <j(A/^ T ) c i a + -^0 an d the partition function is 



Z = J [dc][dJ]exp( -J drL[c,c f ] 



(19) 



Here is the inverse temperature. The integral occur- 
ring in the above expression cannot be calculated in a 
simple way due to the quartic terms in the Lagrangian. 
However, these quartic terms can be eliminated by a 
Hubbard- Stratonovich transformation based on the in- 
tegral identity 



exp(VX 2 ) oc J dYexp(-Y 2 /W + XY). 



(20) 



In our case the quantities corresponding to X are 
E/t c v and E ff c i,A+v Therefore we have to in- 
troduce 2iV bosonic fields: 0i and Xi,i+i- We note that 
the fields <pi and Xi,i+i correspond 2 to site- and bond- 
centered densities, respectively. 

Adding the appropriate terms to the Lagrangian gives 



L[c,c\cf),x} = i[c,c f ] 

„ Ufy. t c . _2 

+ [/^ 2 +l „lz.w 3 

71 ^ — 1 t/ ^ — + 



(21) 



The new Lagrangian is quadratic in the fields, however 
the fermionic fields are coupled to the bosonic ones. The 
explicit form of the Lagrangian is 

i[c,c f ,(/),x] = p>(X +i &) c L* + 7jlXM+i| 2 

+ [(Xi,i+l - t)c\ iCr c i+lt(7 + h.c] + ^0 2 - i^X. (22) 

The partition function is obtained by integrating over c 
and c': 



Z[4>, X ] = J [dc] [dct]exp (-£ At L[c, c\ 0, 



(23) 



Writing it in the form Z[<f>, x] = ex p(~<Scff [0, x]); this 
defines the effective action. The free energy can be ex- 
pressed in a usual way via Z[<j>, x] as 



F[<f>,x\ = -V0 H Z \M)- 



(24) 



As mentioned in the previous subsection, if spatial os- 
cillations occur in the system, they are expected to ap- 
pear with wave number k* = 2irf. Thus we may expect 
spontaneous trimcrization at one-third filling. Therefore 
we suppose that the boson field 4>i takes three differ- 
ent values depending on whether i = 31, i = 31 + 1 or 
i = 31 + 2 with integer I. They will be denoted as 0i, 02, 



and 03 . Similar assumption holds for the fields Xi,i+ij 
too. The three values are Xii X2, and x 3 - The lattice is 
thus decomposed into three sublattices. 

To get a real free energy the fields Q (a — 1, 2, 3) are 
redefined by continuing to the complex plane (i0 Q — * Q ). 
The free energy can then be written as 

n , n , r 1 A f Nn 2 Nn j2 Nn \ 

F({4> a }, {X a }) = 3 X) ( 4j*« + Iff ^ - — 4° ) 

+ ^(i?(fc)-l/3), (25) 



and the summation for k has to be performed over the 
reduced Brillouin zone which is now one-third of the orig- 
inal one (k runs from — it/ 3 to tt/3) and E{k) is the energy 
spectrum of a single fermion coupled to the boson fields. 
It is the eigenvalues of the Hamiltonian 



-i+2 



H = Y, [{(Xi - t)4b i+1 + ( Xa - t)b\ +1 

i 

+ (X3 ~ t)c\ +2 a i+3 + h.c.) 
+ <h.a\ a i + fat>l +1 b i+1 + 3 ct +2 c i+2 



(26) 



where the operators a, b and c belong to different sub- 
lattices. In order to determine the one-particle spectrum 
one has to diagonalize Hamiltonian (|26j) in momentum 
space. Therefore, we are looking for the eigenvalues of 
the matrix 



01 (Xi - t)c 

(Xi - t)e ik 2 

-ik 



(X3-t)e* k 
(X2 - t)e~ lk 



SX3-t>~ lk ( X 2-t)c 



ik 



(27) 



The energy spectrum E(k) has three branches corre- 
sponding to the three (real) solutions of the third-order 
eigenvalue equation: 



# 1 (fc)-(0 1 +0 2 + 3 )/3 = 

- sign(Q)vTP| cos 
E 2 ,i{k) - (0i + 02 + 3 )/3 = 
- sign(Q)yjp| cos 



Q 



(|p|)3/2 



(28a) 



(28b) 



r 1 -if 

- cos 


Q 




(|P|)3/2 


.3 V 





Here P and Q are the parameters in the eigenvalue equa- 
tion when transformed to the form E 3 (k)+PE(k)+Q = 
with E(k) = E(k) - (0i + 02 + 3 )/3. 

The minimum of the free energy of a system with such 
a spectrum cannot be evaluated quite generally. Fortu- 
nately, we are only interested whether the three <p a and 
Xa are different or not. This analysis can be carried out 
easier in terms of the linear combinations: 



A0i 
A0 2 



i-i + 02 + 3 )/3, 

i-i - 2 )/2, 

ii + 02 - 20 3 )/6, 



(29a) 
(29b) 
(29c) 
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and similar definitions for x, Axi and A-X2- One finds 
that although the free energy has an extremum at A0i — 
A02 = and A%i = A%2 = 0, the free energy of the uni- 
form state is not a local minimum. Thus a density-wave 
has to appear in the system. We can conclude that the 
SU(rt) Hubbard model is unstable against the Heisenberg 
coupling in the large-n limit and exhibits inhomogeneous 
spatial ordering for U > 0. 

^From this analysis alone - due to the rather compli- 
cated one-particle spectrum - we cannot decide whether 
the system is dimerized, trimerized, or some other pe- 
riodicity occurs, and whether the density wave is site- 
centered or bond-centered. It is natural to relate the 
nonuniform phase to the fully gaped excitation spectrum. 
Thus when our previous considerations are taken into 
account, trimerized phase is expected in one-third- filled 
models. This will be supported by the numerical calcu- 
lation. We will also see, that the trimerized phase is not 
a special feature of the systems with integer number of 
electrons per site. It occurs in one-third-filled system for 
arbitrary n > 3. 



For convenience the number of sites in the chain was al- 
ways even. Moreover, since we expect dimerized, trimer- 
ized or tetramerized phases depending on the commensu- 
rate filling p/q, the number of sites iV was always taken 
to be an integer multiple of q. 

When a doubling of the lattice periodicity of the 
ground state is indicated by a finite peak in \s(k)\ at 
k* = 7r, a truly dimerized phase gives equal finite values 
for D a and Q s and vanishing T s and P s in the N — > oo 
limit. Stronger and a weaker bonds alternate along the 
chain. When the peak in \s(k)\ appears at k* = 2ir/3 
and a trimerized phase is expected, T s should be finite 
and Q s should vanish. Symmetry considerations imply 
that two equally strong bonds are followed by a weaker 
or stronger bond in this case. In a tetramerized phase, 
the peak in \s(k)\ appears at k* — tt/2, D s , T s , and Q s 
may be finite in the N — > oo limit and only P s vanishes 
necessarily. 



IV. NUMERICAL STUDY OF THE SPATIAL 
INHOMOGENEITY 



A. The numerical procedure 



In this section we present our numerical results ob- 
tained by the DMRG 34 method for the length depen- 
dence of the block entropy sn(1) and its Fourier trans- 
form s(k) to relate them to the number and position of 
soft modes when the model is critical, or to spatial inho- 
mogeneity of the ground state for gapped models. 

The spatial modulation of the ground state can be a 
site- or a bond-centered density wave. A site-centered 
density wave would manifest itself in an oscillation of 
the von Neumann entropy of single sites, Si, with i = 
1, . . . , N or in the local electron density defined by 



(30) 



where |\&gs) is the ground-state wavefunction. The wave 
number of the charge oscillation can again be determined 
from peaks in the Fourier transform of Si or (n*) denoted 
as S\{k) and n{k), respectively. 

The existence of a bond-centered density wave can be 
demonstrated by studying the variation of the bond en- 
ergy or the two-site entropy along the chain. To avoid 
boundary effects we have calculated the difference of two- 
site entropies in the middle of the chain, between first, 
second, third and so on neighbor bonds: 

D S {N) =S N /2,N/2+l ~ SN/2+1.N/2+2 , (31a) 

T S (N) —S N / 2 , N/2+1 — s N/2+2, N/2+3 7 (31b) 

Qs{N) =SjV/2,JV/2+l - s N/2+3,N/2+A i (31c) 

Ps(N) =S N / 2 , N/2+1 - s N/2+4,N/2+5 ■ (31d) 



The numerical calculations presented in this paper 
have been performed on finite chains with open bound- 
ary condition (OBC) using the DMRG technique and the 
dynamic block-state selection (DBSS) approach i 35 ' 36 We 
have set the threshold value of the quantum information 
loss x t° 10~ 5 for n = 3, 4 and to 1CP 4 for n = 5 and 
the minimum number of block states M m ; n to 256. In 
spite of the large number of degrees of freedom per site 
in the n = 5 case the entropy analysis allows one to 
study this problem as well. The ground state has been 
targeted using four to eight DMRG sweeps until the en- 
tropy sum rule has been satisfied. The accuracy of the 
Davidson diagonalization routine has been set to 1CP 7 
and the largest dimension of the superblock Hamiltonian 
was around three millions. As an indication of the com- 
putational resources used in the present work we note 
that the maximum number of block states was around 
1600 for n = 3 and 900 for n = 4 and 5. 

The large- N limit of the entropies and amplitudes of 
the peaks in the Fourier spectrum can be obtained if 
appropriate scaling functions are used. In a critical, gap- 
less model, in leading order, these are expected to scale 
to zero as l/N while in a noncritical model the scaling 
function depends on the boundary condition. Therefore 
for any quantity A the finite-size scaling ansatz 



A(N) = Ao + a/N 13 



(32) 



is used to evaluate the data obtained with OBC, where 
Aq, a, [3 are free parameters to be determined by the fit. 
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B. The numerical results 

1. Models with q — n 

The 1 / n-fillcd case (q = n) has already been considered 
in Ref. [fj and some of the results were listed in Sec. II. 
As has been shown in Fig. [TJ sjv(Z) oscillates with period 
n for finite systems. These oscillations are due to the 
soft modes located at wave numbers k* — 2-K/n. Taking 
every nth value only, sjv(0 can be fitted accurately using 
([3]) for relatively short chains already if U is large. After 
a proper finite-size scaling, the fit gives c = n — 1, as 
expected. 

Taking the Fourier transform of sn(1), besides the large 
positive peak at k = 0, additional negative peaks are 
found at the positions of the soft modes, at k* = it, 2ir/3, 
7r/4, and 2ir/5 for n = 2,3,4,5, respectively. Their am- 
plitude vanishes, however, in the N — > oo limit. As a fur- 
ther check that there are neither site- nor bond-centered 
oscillations in the ground state we have analyzed s,, (rij) 
and Si,i+i. All Fourier components of these quantities 
scale to zero in the thermodynamic limit. As an exam- 
ple, the finite-size dependence of n(k*) for n = 3 and 
n = 4 is shown in Fig. [6] at U = 10. 



.x 10" 



o 


SU(3),f=1/2, k* 


=ji 





SU(3), f=1/3, k 


=2n/3 


□ 


SU(3), f=2/5, k' 


=4n/5 




SU(4), f=1/2, k* 


=71 


+ 


SU(4), f=1/4, k* 


=K/2 


* 


SU(4), f=1/3, k* 


=2n/3 




0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 
1/N 



FIG. 6: Finite-size dependence of n(k*) for various n and 
fillings for U = 10. The solid line is the finite-size-scaling fit. 

When the finite-size scaling of D s , T s , and Q s is ana- 
lyzed one finds that they all vanish in the thermodynamic 
limit as shown in Figs. H and [9] for U = 10. All this 
shows that the ground state of the 1/n-filled SU(rt) Hub- 
bard model is spatially homogeneous, the translational 
symmetry is not broken. 



2. Models with q > n 

We have chosen as an example n — 3 and / = 2/5. As 
seen in Fig.[3^a), the block entropy sjv(Z) oscillates with 
period 5. When every fifth data points are fitted to |3]), 
c = 3 is obtained. This indicates that the model remains 
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2- 
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SU(3),f=1/3 
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* 


SU(4),f=1/3 


+ 


SU(4),f=1/4 
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FIG. 7: Finite-size dependence of D s for various n and fillings 
for U = 10. The solid line is the finite-size-scaling fit. 
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FIG. 8: Finite-size dependence of T 3 for various n and fillings 
for U = 10. The solid line is the finite-size-scaling fit. 



critical for finite U as well. The finite peak in \s(k)\ at 
k* = 4ir/5 is due to soft modes. The amplitude of the 
peak disappears in the N — > oo limit, the ground state 
of the system is uniform. This is confirmed by the calcu- 
lation of s(k*), n(k*), D s , T s , and Q s shown in Figs. 
El El [H and [9l respectively. There is neither a site- nor 
a bond-centered oscillation in the occupation number or 
bond strength. 



3. Models with q < n 

One realization of this condition, the half-filled case 
for n > 2 has been studied by us earlier^ It was found, 
as shown in Fig. [U that the block entropy oscillates with 
period 2 for any n > 2. The peak in \s(k) \ at k* = ir does 
not vanish in the thermodynamic limit (see Fig. [5]). In 
agreement with this D s and Q s are finite and converge 
to the same value as shown in Figs. [7] and [£l while T s 
vanishes (see Fig. HJ. The same behavior is found in the 
n = 5 model at half filling, as seen in the upper panel of 
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FIG. 9: Finite-size dependence of Q s for various n and fillings 
for U = 10. The solid line is the finite-size-scaling fit. 
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FIG. 10: Finite-size dependence of D a ,T a , Q a , and P s for 
the half-, one-third-, and quarter-filled one-dimensional SU(5) 
Hubbard model at U = 10. 

On the other hand, s\(k*) and n(k*) vanish in the 
thermodynamic limit (see Fig.|BJ). The translational sym- 
metry of the Hamiltonian is broken, the ground state of 
the half-filled model is dimerized. Stronger and weaker 
bonds alternate along the chain. 

In the one-third-filled case of the n — 4 and n = 5 
models (q < n) sjv(i) oscillates with period 3 and the am- 
plitude of the Fourier component s(k* = 27r/3) remains 
finite even for N — > oo (see Figs. [2] and 2]). This indicates 
that the spatial periodicity is tripled in the ground state. 
This is corroborated by our results shown in Figs.[7l[8j[9j 
and the middle panel of [TUl D s and T s scale to the same 
finite value while Q s , s\(k*), and n(k*) vanish. In the 
ground state, two bonds of equal strength are followed 



by a weaker or stronger bond. 

As a last example we have studied the quarter-filled 
SU(5) Hubbard chain. We found that |s(fc)| scales to a 
finite value at k* — it /A as shown in Fig. 2J All Fourier 
components of the site entropy and local charge density 
vanish for long chains, while D s , T s , and Q s scale to fi- 
nite values. Only P s scales to zero, as shown in the lower 
panel of Fig. [TUl The ground state is a bond-ordered 
tetramerized state. Fig. [TTJ shows schematically the pe- 
riodic modulation of the bond strength along the chain 
for half-, one-third, and quarter-filled models. 



Half-filled 



One-third-filled 



Quarter-filled 



Dimerized 

JO 
Trimerized 

X> 

Tetramerized 



FIG. 11: Schematic plot of the local bond strength in half-, 
one-third, and quarter-filled SU(5) Hubbard chains. 



V. CONCLUSION 

To study the role of multiparticle umklapp processes, 
we have treated the one-dimensional SU(n) Hubbard 
model analytically by bosonization approach and numer- 
ically using the DMRG method for n = 3, 4, and 5 for 
commensurate fillings / = p/q where p and q are relative 
primes. 

Our results confirm that umklapp processes play es- 
sentially different role depending on the relationship be- 
tween q and n. When q = n (this is the case in the 
1/n-fflled case) the charge and spin modes are not cou- 
pled, the umklapp processes open gap only in the spec- 
trum of charge modes. The system remains critical with 
ri — 1 gapless spin modes, the central charge is c = n — 1, 
and the translational symmetry of the Hamiltonian is not 
broken in the ground state. 

When q > n, the leading-order umklapp processes are 
forbidden in the model with local interaction by Pauli's 
exclusion principle. The model is equivalent to an n- 
component Luttinger liquid with c = n and the ground 
state is spatially uniform for U > 0. 

When, however, q < n the charge and spin modes are 
coupled by the umklapp processes and gap opens in the 
spectrum of all modes. Even more interestingly, a spa- 
tially nonuniform ground state emerges whose periodic- 
ity depend on the filling. Half-filled models develop a 
dimerized ground state, trimerized state appears in one- 
third-filled models, the ground state is tetramerized in 
quarter-filled models. Other periodicities would proba- 
bly be found at other fillings. Our findings are summa- 
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rized in Table HI1 which can be compared to the analytical 
results given in Table HI 



TL 


p/q 


c 


periodicity 


h 


q — TL z 


1 /9 
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q < n 3 


1/2 




dimerized 


TV 
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dimerized 


TT 
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trimerized 


2tt/3 
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1/2 




dimerized 


7T 


5 


1/3 




trimerized 


2tt/3 


5 


1/4 




tetramerized 


tt/2 


q > n 3 


2/5 


3 


uniform 


4tt/5 



TABLE II: Central charge and spatial inhomogeneity for the 
p/q- filled SU(n) Hubbard chain, k* in the last column gives 
the wave number of soft modes when the model is critical 
while it gives the wave number of the nonuniform ground 
state when the model is gapped. 



We emphasize that our calculations were performed at 
a relatively large value of U, where the nonuniformity of 
the ground state is well developed, and the finite value 
of the dimer, trimer or tetramer order parameter can 
easily be detected. We conjecture, based on our earlier 
calculations^ that the critical value U c above which the 
nonuniform phase appears, is U c = 0- 
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